#include "ELAPLA.h"
#include <string.h>

void multi(double *a, double *b, double *c)
{
    double temp = 0.0;
    int N = MATRIX_SIZE;
    memset(c, 0, sizeof(double) * MATRIX_ALLOCATE_SIZE);
    for (int i = 0; i < MATRIX_ITER; i++)
    {
        for (int j = 0; j < MATRIX_ITER; j++)
        {
            for (int k = 0; k < MATRIX_ITER; k++)
            {
                c[i * N + j] += a[i * N + k] * b[k * N + j];
            }
        }
    }
}
